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We study the quantum phase transitions in two-dimensional arrays of Josephson-couples junctions 
with short range Josephson couplings (given by the Josephson energy Ej) and the charging energy 
Ec- We map the problem onto the solvable quantum generalization of the spherical model that 
improves over the mean-field theory method. The arrays are placed on the top of a two-dimensional 
electron gas separated by an insulator. We include effects of the local dissipation in the presence of 
an external magnetic fiux / = "l>/"3>o in square lattice for several rational fluxes / — 0,^,^,j and 
|. We also have examined the T = superconducting-insulator phase boundary as function of a 
dissipation qq for two different geometry of the lattice: square and triangular. We have found critical 
value of the dissipation parameter independent on geometry of the lattice and presence magnetic 
field. 
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; I. INTRODUCTION 

Over the past several years Josephson junction arrays 
(JJA) have gained increased interest since the display 
quantum mechanics on macroscopic scale. In the last 

[ years, due to development of the micro-fabrication tech- 
niques, it became possible to fabricate Josephson junc- 

' tion arrays with tailor specific properties. In these sys- 
tems the superconductor-insulator (SI) transition can be 
driven by quantum fluctuations when the charging energy 
Ec becomes comparable to the Josephson coupling en- 
ergy ii;j i^i^ii?i"^ii?ifii7iff It has been established that in arrays 
which are in the superconducting state, but with value 
Ec/Ej close to critical value, a magnetic field can be 

' used to drive array into the insulating state.™ There are 
experimental possibilities to fabricate different structures 

I of the JJA like squar o^ft'^^'^^ and triangular- 4 networks 
which show that the critical value Ej j Ec can change its 
value depending on the geometry of the lattice. Zant et 
al. showed experimentally that not only geometry can 
influence on value Ec / E,j in JJA but also external mag- 

\ netic field can lead system to the phase transitioniiSii^ In 
a JJA an applied gate voltage relative to the ground plane 

I Vg introduces a charge frustration^ The combination of 
charge frustration and Coulomb interaction results in the 
appearance of various Mott insulating phases separated 
by regions of phase coherent superconducting state. 

It has been understood early that dissipation can 
be capable of driving an SI transitions. Problem 
how dissipation can be described on quantum me- 
chanical level was firstly addressed by Caldeira and 
Leggettfl^ In this formalism dissipation introduces damp- 
ing of phase fiuctuations that is inversely proportional 
to the resistance of the environment. Quantum phase 
model of JJA was formulated in terms of an effective 



actioufiS^ Various theoretical methods have been ap- 
plied in effects of Ohmic dissipationiS'ifi'SfliSl such as 
coarse graining jSSi^ variationaU^'^^'^^i^^'^^^ and renor- 
malization group approaches piSiSiiS2i21 The dissipation 
due to quasi-particle tunneling in JJA2£ was investigated 
by means of the mean field calculations variational 
approaches'^ and Monte Carlo simulations Phase tran- 
sitions of dissipative JJA's in magnetic field were ana- 
lyzed by Kampf and Schon^ and relied upon mean-field 
approximation mapped into tight-binding Schrodinger 
equation for Bloch electrons in magnetic field. It has 
been shown that at a fixed temperature we can observe 
a phase transition if we vary the magnetic field. The 
similar problem was solved by Kim and Choi^^ based on 
effective Hamiltonian obtained by the Feynman path in- 
tegral formalism. Authors claim that especially in low 
temperatures variational method is not precise enough 
to perceive such a subtle transition. 

Dissipation may arise from coupling with the substrate 
by means of the local damping modeli^ These studies in- 
deed revealed the existence of a critical value of the sheet 
resistance which seems to agree with experimental results 
in JJA24 and thin filmsj24 On the other hand some ex- 
perimental studies discloses that the values of the critical 
resistance can show wide variations quite the contrary to 
the predicted universal value close to /i/4e^i^ Further- 
more, previous theoretical calculations rehed upon varia- 
tional and mean field approximations, which usually are 
not expected to be reliable at T = and be capable to 
handle spatial and quantum fiuctuations effects properly, 
especially in two dimensions. 

It has been shown^S that there is a possibility of ex- 
istence four phases in superconducting junction arrays 
with dissipation and the phase diagram depends on ratio 
Ej /Ec and dissipation parameter ao . One is insulating, 
when both Cooper pairs and single electrons are frozen 
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{Ej/Ec and ao small), mixed, when superconducting 
long-range coherence and single electrons tunnehng co- 
exist {Ej I Ec and ao large) , and for intermediate values 
of the parameters we can obtain {Ej / Ec small and ao 
large) Cooper pairs remain frozen, but single electrons 
are free. In the opposite limit {Ej/Ec large and ao 
small) we have superconducting long-range phase order 
and there is not any dynamics of the single electrons. 

Fact, that dissipation could play important role in solid 
state physics appeared recently from the high-Tc point 
of view. Relation obtained by Homes et aWM proves that 
the characteristic timescale for dissipation could not be 
shorter than in the high-Tc superconductors. Relaxation 
time Tdiss ~ h/ksTc (Planckian dissipation) is an essence 
of the Home's law and there is evidence that, quantum- 
critical nature of the system could be present even in the 
normal state. It means that below this timescale we have 
purely quantum mechanical behaviour and energy could 
not be turned out into heat. Conductivity in the normal 
state is tied with the relaxation time relation (Tnormai ~ 
Tdiss and it indicates that (Jnormai should be as small as 
it is allowed by Planck's constant. 

Realization of a quantum computer crucially depends 
on our ability to create and hold coherent superposition 
states, so decoherence presents the most fundamental 
trammel. Especially coupling between different devices 
and environment achieves to dissipation, and hence deco- 
herence which leads to exponential decay of superposition 
states into incoherent mixturesi^^ Both in superconduct- 
ing qubits, based on superconducting interference devices 
and in single-pair tunneling devices the Josephson junc- 
tion is a key element and it is the dissipation of the junc- 
tion that sets the limit on decoherence time4fi 

To purpose of this paper is to study local dissipation 
effects in JJA in an analytical way to refine the phase 
diagram of the system for different geometry of the lat- 
tices and in presence of the external magnetic field. We 
consider a network of the Josephson junction arrays with 
tunable dissipation which is placed on the top of a two- 
dimensional electron gas (2DEG) separated by an insu- 
lator. We drop nonlocal charging and dissipative terms. 
The problem we would like to address is then: What is 
the effect of the competition between magnetic, geomet- 
ric and quantum effects on the ground state ordering in 
the two-dimensional Josephson arrays? To analyze 2D 
JJA beyond mean field approximation we employ the 
path-integral formulation of quantum mechanics and a 
quantum spherical model approach accurately tailored 
for the JJA systems. This formalism allows then for ex- 
plicit implementation of the local dissipation effects and 
magnetic field into ours considerations. Most theoretical 
studies investigated the simple square lattice geometry of 
JJA. Another structures were analyzed by Monte Carlo 
simulations and mean field calculations in magnetic field 
in the context of phase transitions^. On the other hand 
Granato and Kosterlitz claim that transition in 2D array 
with differential geometry can be described by classical 
Ginzburg-Landau- Wilson effective free energy with two 



complex field4^ We analyze the quantitative changes in 
the phase diagram due to two different geometrical JJA 
structures without external magnetic field. We consider 
infiuence of the magnetic field on square lattice because 
many different properties of an array depend on fiux pa- 
rameter / — p/qM^ 

The paper is organized as follows. In the next sec- 
tion we introduce the model. In Sec. lie we formulate 
the problem in terms of the effective action of quantum 
spherical model. The various zero temperature phase di- 
agrams are then studied in Sec III. Finally, in Sec IV we 
summarize. 



II. MODEL 

We consider a two-dimensional Josephson junction ar- 
ray with lattice sites i, characterized by superconducting 
phase (pi in dissipative environment. Possible experimen- 
tal realization of the 2D JJA is shown on the FigQ] The 
array can be formed by thin square superconducting is- 
lands of size L. The separation d between islands must 
be small enough to allow for the Josephson interactions. 
The variable dissipation is introduced by coupling with 
two-dimensional electron gas (Ohmic bath) which is lo- 
calized within distance s. The corresponding Euclidean 
action in the Matsubara imaginary time r formulation 
(0 < T < l/k-eT = P), with T being temperature and /cb 
the Boltzmann constant (ft — 1) 



where 
Sc 



S = Sc + S:i + Su, 



V / drJy {1- cos [(/)i(T)-(/)j(r)]}, 

^ E / d^dr'a (r - r') (r) - 0, (r')]' .(2) 



(1) 



The first part of the action Eq. ^ defines the electro- 
static energy, where self charging energy parameter is 



Ec 



2Co' 



(3) 



The second term is the Josephson energy Ej {Jij = Ej 
for \i — j\ = \d\ and zero otherwise). Third part of the 
action Sn describes the local dissipation where a (r — r') 
is dissipative kernel. For a local dissipation effects Fourier 
transform (with respect to imaginary times) of the kernel 
Eq. © isiiSa 



a (Wn) = — , 



(4) 
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from the microscopic model Eq. Several studies of 

JJA have followed this way, also known as coarse grained 
approach first developed by Doniachii Most of existing 
analytical works on quantum JJA have employed dif- 
ferent kinds of mean-field-like approximations which are 
not reliable for treatment spatial and temporal quantum 
phase fluctuations which quantum spherical approach 
captures. Formulation of the problem in terms of the 
spherical model^^*^ which was initiated by Kopec and 
Jose^ leads us to introduce the auxiliary complex fleld 
tpi whose expectation value is proportional to the Si {(p) 
vector defined by 



Figure 1: Schematic view of the 2D JJA in dissipative envi- 
ronment . 



where dimensionless parameter: 



Rq 

Rq ' 



1 



(5) 



describes strength of the Ohmic dissipation. Here the 
i?o is the shunt to the ground and is interpreted as the 
shunt resistance present in many experiments. As we can 
see the dissipation part of the action Eq. (EJ breaks the 
27r-periodicity in the phase variables since it allows for 
continuous charge fluctuations. In proximity-coupled ar- 
rays, dissipation can correlate the phase of a single island 
in different time. 



A. Effect of the applied magnetic field 

The perpendicular magnetic fleld B given by vector 
potential A enters action Eq. (EJ through a Peierls factor 
according to 



J.,exp( — 



A-d\ 



(6) 



where $o = 27rc/e is a quantum of magnetic flux piercing 
a 2D lattice plaquette and integral runs between center of 
grains and rj. Thus, the phase shift on each junction 
is determined by the vector potential A and in typical 
experimental situation it can be entirely ascribed to the 
external field B. The presence of B induces vortices in 
the system described by the fiux parameter /(= <i>/<i>o — 
Ba^/^o where a is area of an elementary plaquette). Our 
special interest are the values of the magnetic field when 
fiux parameter f — p/q represents rational values. 



B. Method 

To study the JJA model we can use a description in 
terms of an effective Ginzburg-Landau functional derived 



S. (0) = [S^ (4>) , 5? m ^ [cos {ct>^) , Sin (00] • (7) 

Natural consequence definition (0) is the following 
rigid constraint 

|S, (0)|2 = [St {4>)f+[Sl {4>)f = cos^ (0,)+sin2 {<t>i) ^ 1. 

(8) 

The relation in Eq. l|Hl implies that a weaker condition 
also applies, namely: 



(9) 



Using the Fadeev-Popov method with the Dirac 5 func- 
tional we obtain: 



Z = 



e-5c+°M Yl S [3?V, (r) - (0 (r))] 

i 

xS[^^,ir)-Sf{^{r))]. (10) 



It is convenient to employ the functional Fourier rep- 
resentation of the S functional to enforce the spherical 
constraint Eq. ©: 



S[xiT)] 











J —ioo 


2TTi 


exp 





drX (t) X (r) 



(11) 



which introduces the Lagrange multiplier A (r) thus 
adding an quadratic term (in field) to the action Eq. 
^ . To evaluation of the effective action in terms of the 
-0 to second order in -01 (t) gives the partition function 
of the quantum spherical model 



--QSM 



S \^ (^)I' - e-'5QSAW (12) 
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where the action of the eflective nonlinear a model reads: 



5qsa[^,A] = / drdr' {[J,,{t)5{t-t') 

+ yV-j^ (r, r') - A (r) 5,,5 (r - r')] 

+ N\{t)6{t-t')}. (13) 



Furthermore 



2o 



(15) 



where action 5c+d ['^] is just a sum electrostatic and dis- 
sipative term in Eq. After introducing the Fourier 

transform of the field 



^0 



= >V(t,t')%, 



(r')] -5c+D [ 



(14) 



is the phase-phase correlation function with statistical 
sum 



<Pi ('^) = X! X! "^^P [~ (^^"■^ ~ ^''»)] (-^^^ 

k n— — oo 



with a;„ = 27rn//3, (n = 0, ±1, ±2, ...) being the Bose 
Matsubara frequencies, we can write expression Eq. Ijl4|l 
in form 



+ 00 r 



exp ■ 



— y 

MR ^ 



k 

+ 00 



n— — OO 



1 2 "0 I 



»k,«(/'k,i 



''k.ri 



(17) 



Using Hubbard-Stratonovich transformation the phase- 
phase correlation function reads: 



mr,.^^^{ ^^':X^}' Y (18, 

The low temperature properties of the expression 
(ujn) lead to critical value — 2 [see Appendix 
A]. Finahy, for small frequencies and ao < 2 inverse of 
the correlation function Eq. itTsll becomes: 



for cj„ 7^ and (w„) = otherwise. In the thermo- 
dynamic limit (iV oo) the steepest descent methods 
becomes exact; the condition that the integrand in Eq. 
II12II has a saddle point A (t) — Aq, leads to an implicit 
equation for Aq: 



^ k n^^O 



where 



G ^ (k, w„) = Ao - Jk 



2-K 



(20) 



(21) 



with Jk as Fourier transform of the Josephson couplings 



Jij. 



III. PHASE DIAGRAMS 

A Fourier transform of Eq. Ijl3ll in momentum and 
frequency space enables one to write the spherical con- 
straint Eq. ll2?11l explicitly as 



'4 



4-00 



SEc ^" ' 



2ir 



0J„ 



(22) 



As is usual in a spherical model, the critical behavior and 
the phase transitions boundary crucially depends on the 
spectrum given by density of states p(^) and is deter- 
mined by the denominator of the summand in the spher- 
ical constraint Eq. II22II . To proceed, it is desirable to 
introduce density of states for two dimensional lattice in 
form 



(23) 



where Jk is energy dispersion. Problem of computing 
density of states p^^^ (^) for two-dimensional square lat- 
tice with uniform magnetic flux per unit plaquette re- 
duces to the solution of Harper's equation^, e.g. rel- 
evant for tight-binding electrons on 2D lattice^. An- 
alytical results for density of states for square lattice 
were presented recently®, and based on analytically solv- 
ing Harper's equation and receiving dispersion relation 
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J^p/q- Influence of the local dissipation effects will be 
considered on triangular lattice without magnetic field. 
In that case we can use definition Eq. ^2?>\ but the de- 
pendence on the wave vector is different and could be 
described by expression 



cos [kx) + 2 cos 



2 



cos 




(24) 



Closed formulas for the density of states are placed in 
Appendix B. The systems displays a critical point at Aq = 
(k = 0, w„ = 0) = max [J (k = 0)]. This fixes the 
saddle point value of the Lagrange multiplier, A sticks to 
that value at criticality (A = Aq) and stays constant in the 
whole low temperature phase. By substituting the value 
of Ao to Eq. 1I22I1 and after performing the summation 
over Matsubara frequencies, in T ~> limit we obtain 
the following result: 



1 = 



+ 00 



p{0 



m - 



2Ec 



X In 



2lT 



.2Tr ) 



2Ec 



aq 
2tt 



aq \ 
.2itJ 



Jq-^Ej 
2Ec 



(25) 



It is easy to see that by specifying density of states p (^) 
Eq. lf23jl with the Coulomb energy Ec Eq. ((HJ for given 
superconducting network without and in presence of the 
external magnetic field, we are able to study the zero tem- 
perature JJA phase diagram. The solutions and bound- 
ary cases above equation we will examine in the next 
subsections. 



A. Small ao limit 

In the limit ao ^ expression II25|I reduces to zero- 
temperature critical line in absence dissipation effects: 



1 



2Er 



Jo - ^Ej ' 



(26) 



which is according with previous calculations^. When ao 
is nonzero, but still small, the situation changes due to 
a dissipation is a factor which produces additional quan- 
tum fiuctuations. In consequence for small ao the critical 
value Ej/Ec decreases as: 



Ej/Ec 



42 
^0 



2 



2 Ai 
yA'^ 
2 Ai 

where corresponding coefficients reads: 



2 Ai 



^0 



a^ + O (a^) 



Aq = V2 



-f 00 



^Jo/Ej - e 



(27) 



(28) 
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Figure 2: a) Zero temperature phase diagram for the total 
charging energy Ej/Ec vs parameter of dissipation qq for 
square lattice. Insulating state is below the curve. From 
the top we have p/q = \ and b) Nonmonotoni- 

cal dependence the critical value of the inverse Coulomb en- 
ergy Ej /Eq^^ for several ratio of the magnetic fluxes p/q ~ 

0, 21 3 ' 4 ^^'^ 6 • 



DOS 






„u 

Pl/2 


Pl/3 


„u 

Pl/4 


Pl/6 


Ao 
Ai 


1.01087 
1.92619 


1.28576 
5.82281 


1.34085 
4.83852 


I. 65397 

II. 8065 


2.06193 
19.1268 


3.78672 
175.717 



Table II: Factors of the expansion critical values Ej /Ec for 
small dissipation parameter qq (Eq. 1281 129|l . 



Ax = 



Jo/Ej - C 



(29) 



The numerical values of the factors ^0 and Ai are clas- 
sified in the Table For small ratio Ej/Ec there 
is no chance to mobility both Cooper's pairs and sin- 
gle electrons. If we increase the Josephson energy (or 
decrease Coulomb energy) we induce phase transitions 



6 



between insulating and superconducting phase. We have 
global coherence state but in which there is no single 
electrons dynamics. The critical values of the inverse 
Coulomb energy Ej/Eq"^* for ao = are depicted in 
FigEb) for several values of the magnetic field and are 
simply equal the first coefficient of the expansion Eq. 
ll27jl : Ej/E^'^^l^^^^^ = ^0- We note the non-monotonic 
dependence of the Coulomb energy on the magnetic fiux 
ratio p/q. 

B. Critical value ao 

At zero temperature dissipation suppresses quantum 
fiuctuations entirely and drives system to a global co- 
herent state. Due to the fact, that correlation function 

Eq. II18II for low temperatures: (ijJ„) ~ |wn|°" ^ be- 
comes divergent for ao > 2 [see Appendix A], the critical 
lines are cut at this point what is depicted on phase dia- 
grams in FiglSK and FigEl This boundary value ao does 
not depend on magnetic field and on the geometry of the 
network. It seems to be universal for different types of 
lattices without and in presence of the external magnetic 
field. The system behaves as if it were classical because 
of the big contribution to the action Eq. (01 which is 
generated by large values of the dissipation parameter 
ao. Dissipation can correlate the phase on a single is- 
land in different time and this correlation has the biggest 
impact on global coherent state when ao is greater then 
critical value. Similar behaviour of the phase diagrams 
with critical value of the dissipation parameter was pre- 
dicted by several authorsii^iSiii^ Theory developed by 
Chakravaraty et alm^ reveals, that phase transition takes 
place at point acrit = l/d where d is the dimension of 
the system. The insulating phase disappears for ao > 2 
in model considered by Wagenblast et alM-, the authors 
claim that dissipative processes completely suppress the 
phase fiuctuations. Kampf and Schon^ using variational 
procedure showed that different mechanisms: Ohmic and 
quasi-particle damping lead to different critical values of 
ao. In magnetic field the dependence of critical temper- 
ature on several ratio of the magnetic fiuxes is periodic 
with period 1 and symmetric around p/q = 1/2^ Note 
that phase diagrams for different lattice geometry and 
in presence of a magnetic field with effects of the local 
dissipation has not been presented. 

IV. SUMMARY 

We have studied the ground phase diagram two- 
dimensional Josephson junction arrays in quantum 
regime by analytical methods using the quantum- 
spherical approach with exactly evaluated density of 
states for square and triangular lattice. For square lat- 
tice we could take into consideration perpendicular mag- 
netic field which was described by rational magnetic fiux 
f = p/q for a number of values 1/q. Zero temperatures 





1 1 










1.5 












- 








1.2 






^"-^ 1 — 1 




U.y 


- 










0.6 












0.3 






A 








, 1 




1 , 


, 1 





0.3 0.6 0.9 1.2 1.5 1.8 a^2.1 



Figure 3: Zero temperature phase diagram for the normalized 
Josephson energy Ej/Ec vs parameter of dissipation ao for 
square and triangular lattice. Insulating state is below the 
curves. 



phase diagram indicates that for small values ao quan- 
tum fiuctuations destroy the long-range phase coherence. 
The arrays can be in two phases: Mott-insulator phase 
and global coherent state. We can travel between phases 
vary Coulomb energy or dissipation parameter. When ao 
is greater than critical value, the dissipation suppresses 
quantum fiuctuations and the array orders even for small 
ratio Ej/Ec- Why dissipation can restore global coher- 
ent state? If we imagine that we allow electrons tunnel 
between superconducting arrays and metallic substrate, 
remaining Cooper's pairs will become more mobile too. 
Therefore we will not be equal to set a number pair on 
the island. Because superconducting phase 4> and the 
number of Cooper's pair fi follow uncertainty relation 
[4>, n] = 2i then if we cannot say anything about quan- 
tity of Cooper's pair on the island the phase is determined 
and in consequence we obtain global coherent state. In 
case when system is in presence of the magnetic field 
damping is stronger than in the absence. If we vary fiux 
parameter f = p/q we could drive array into insulating 
state. Magnetic field could affect the dissipation because 
it infiuences the resistance of at least the coherent parts 
and thus changes the effective shunting. For small values 
dissipation parameter ao Josephson energy in triangular 
lattice is damped less then for square lattice even in pres- 
ence magnetic field. We can notice that if the quantum 
fiuctuations of the phase superconducting order param- 
eter are important, the dissipation plays decisive role in 
constitute the onset of global superconductivity. 
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Appendix A: LOW TEMPERATURE 
PROPERTIES OF THE CORELLATOR 



can write 



We write correlation function Eq. Ijl8|l in form: 



>V(r) = exp<j--5]3 



1 — cos {u>nT) 



SEc "^n 271- 



(Al) 



and observe that the sum over w„ is symmetric when we 
change uJn —^n- In that case (getting rid of abs) we 



W(t) =exp^ -- 



1 — cos {uJn''') 



(A2) 



noticing then sums over n > 1 and n <\ are the same. 
Now we are sphtting above expression for two parts and 
neglecting the in second contribution (in low temper- 
atures more important is part which contains a;„) we get: 



W (t) = exp 

Evaluating the sums we can write: 

W (t) = exp 



2 ^ 1 

a ^ 1 -.2 1 "0,. 



exp 



-y 



cos [uJuT) 



n>l 27r 



exp < — log 

"0 



2 sin I - in 



and using asymptotic relation (valid for (3 oo) of the harmonic number function H (n) = X]r=i 



1. „/l 



H(a/3) = lna/3 + 7+-a/3-- — +_ — + ...O { -r 



2 ' 12 \al3 J 120 Va/3 
where 7 is Euler's constant, finally we have found 

After Fourier transform correlator becomes: 



(c^„) = exp (^^^ r 



1 I sin I — ] |w„| °fi ^ 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



r 



besides the Euler gamma function F (z) is defined by the 
integral 4i 

poo 

F(z)= / t'-^e-^dt for > 0. (A8) 

and can be viewed as a generalization of the factorial 
function, valid for complex arguments. Finally we can 
see that correlator 

W-^ M ^ \LUn\^-^ (A9) 

at zero temperature diverges for ao > 2. 



Appendix B: DENSITY OF STATES 

In the case of zero magnetic field the density of states 
for the square two-dimensional lattice reads simply 



Po iO 



where 



rK 



K(x) 



tt/2 



e 1 



d<j) 



vi^ 



a;^ sm (h 



(Bl) 



(B2) 



is the elliptic integral of the first kind^ and the unit step 
function is defined by: 



I for x < 



(B3) 



The procedure obtaining closed formulas for density of 
states in presence of the magnetic field based on analyt- 
ically solving Harper's equation. The Harper's equation 
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was studied intensiveljiS2iS2i^ but only on numerical way reads 
and was missed analytic closed formulas for density of 
states in presence of the magnetic field. Only expression 
for case p/q — 1/2 was received by Tan and Thouless and 
also was formulated in terms of the elliptic integralsi^ 
Therefore the density of states for square lattice with the 
magnetic quantum flux per plaquette for value p/q = ^ 



Pi/2 (4) ^ -irPo 



e-4 



(B4) 



It is only one gap- less case instead of p/q = 0. 



Obtaining closed formulas for the next case p/q= 1/3 was more difHcult, and we have got following expression 



1 



4\/2 

X 



{|sec (a + 1) I [e + 1 + Vs) - e (6 - e') - e - 1 - \/3' 
+ sec (a + 1) [e + Ve) - e (e + 2) + e (0 - e + 1 - Vs' 
+ sec (a + [e - 1 + Vs) - e (c) + e - 2) e (e - Ve)] } 



where 



1 

a = - arctan 

O 



For p/q — 1/4 the expression for density of states is straightly given by: 



(B5) 



(B6) 



P?/4 (0 



4+1^2-41 e(8-C^) -e(4 + 2\/2-C 



V4-|?2_4|e (4 _ 2\/2 - e) } 



(B7) 



Finally the most complicated case p/q = 1/6 with six symmetric singularities divide symmetrically on the positive 
and negative part of the axis ^ : 



p?/6io = j^\e-8e+8\v¥^w^p^ 



ie - 6) 



X {|sec (a + I) I [e (e + 1 + \/3) - e (6 - _ _ 1 _ 

+ sec (a + J) [e + \/6) - e + 2) + e (0 - e + 1 - \/3 

e - 1 + x/s) - e (0 + e - 2) e (e - x/e)] } 



where 



+ sec I a H 

6 



a = — arctan 

o 



- Se + 8) + 8^2 _ ^4 



^6 - 12^4 _^ 24^2 



32 



(B8) 



(B9) 



The density of states for triangular lattice with six nearest neighbors we can write in form: 



7r2^^ 



K( 

Kq 



e(c-3) 



where 



Kq 



(3+2^3 + 2^-^ 



0(e+2 ) -e(e + i) 



A^/JT2i[e (e + i)-e(^-3)], 



(BIO) 



(Bll) 



-e(c + i) 



(3 + 2^/3 + 21 - [e + 1) - e - 3)] . (B12) 
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